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Abstract 

The phase structure of four-dimensional simplicial quantum gravity 
coupled to U{1) gauge fields has been studied using Monte-Carlo simu- 
lations. The smooth phase is found in the intermediate region between 
the crumpled phase and the branched polymer phase. This new phase 
has a negative string susceptibility exponent, even if the number of vector 
fields {Ny) is 1. The phase transition between the crumpled phase and 
the smooth phase has been studied by a finite size scaling method. From 
the numerical results, we expect that this model (coupled to one gauge 
field) has a higher order phase transition than first order, which means the 
possibility to take the continuum limit at the critical point. Furthermore, 
we consider a modification of the balls-in-boxes model for a clear under- 
standing of the relation between the numerical results and the analytical 
one. 
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1 Introduction 



To formulate a theory of quantum gravity in four dimensions, many approaches 
have been tried. One of them is a numerical approach with the method of dy- 
namical triangulation. In two dimensions, the quantum gravity can be quantized 
for a central charge c < 1. The method of dynamical triangulation has generally 
been considered to be a correct discretized model, and has given consistent results 
with the analytical approach: for example, the MINBU analysis and the loop 
length distributions 0,^]. Two-dimensional quantum gravity has generally been 
regarded as being a toy model of four-dimensional quantum gravity. Recently, 
a numerical approach with the method of dynamical triangulation in analogy 
with a two-dimensional model has been studied. In the four-dimensional case, it 
has been expected that the phase transition point between the strong coupling 
phase and the weak coupling phase is statistically the second phase transition 
point. Moreover, this point is recognized as the ultraviolet (UV) fixed point of 
the quantum theory of gravity. 

From numerical results, in four-dimensional pure gravity, it is known that 
there are two distinct phases p|-p!0|. For small values of the bare gravitational 



coupling constant, the system is in the so-called elongated phase, which has the 
characteristics of a branched polymer phase. For large values of the bare gravita- 
tional coupling constant, it is in the so-called crumpled phase. Numerically, the 
phase transition between the two phases has actually been shown to be of first 



order 12]. Therefore, it is difficult to construct a continuum theory. In other 



words. It is difficult to define the quantum theory of gravity on a four- dimensional 
triangle (4-simplex) as a simple application of the two-dimensional lattice model. 

Our next step is to investigate the possibilities to extend the model of four- 
dimensional quantum gravity. We have three motivations: (1) a modified model 

in three-dimensional case, which suppresses the vertex order concentration 
(VOC) as the singular sub-simplex, and which can be changed of the phase 
structures, (2) the property of quantum field theory with background metric 
independence [|T^T3], that the manifold may be made stable by adding matter 
fields and (3) the balls-in-boxes model, that gives a scenario for a phase transition 
in four-dimensional simplicial quantum gravity [|l6i[T7[] , and that shows possibil- 
ities to change the phase structure with some modifications of the model [0. 
For the possibility of a continuum theory, at least, we have to find a new phase 
structure that has a second order phase transition point. One modification is to 
introduce gauge matter fields . 

Recently, the phase structure with vector fields has been studied numerically 
[^,^. In the case of a model with vector fields, the phase structure is changed 
drastically and the intermediate phase, the so-called smooth phase, has been 
observed between the crumpled phase and the elongated phase. In this region, 
the string susceptibility exponents (7^4) have negative values. They show that this 
region has a fractal property, and may be smooth compared with the branched 
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polymer region. We thus expect the possibihty of a continuum hmit at the 
critical point between the crumpled phase and the smooth phase. In order to 
investigate the nature of the phase transition, we measure the critical exponent as 
the finite size scaling, and also study the scaling property of the mother boundary 
in analogy with the the two-dimensional case, and we expect that the scaling 
structure also appears in the boundary in the four-dimensional case. 

This paper is organized as follows. In Section 2, we discuss the model of four- 
dimensional dynamical triangulation with some vector fields. In Section 3, we 
show our numerical results concerning measurements of the string susceptibility 
exponents {•jst) and a schematic phase diagram. We thus discuss the statistical 
property of the phase transition between the crumpled phase and the smooth 
phase in the case of four-dimensional simplicial quantum gravity coupled to one 
gauge field {Ny = 1). Furthermore, in section 4, we discuss the scaling property 
near to the critical point. In section 5, we discuss a scenario for the phase 
structure and the phase transition in four-dimensional simplicial quantum gravity 
coupled to matter fields. Finally, we discuss the possibility of a continuum limit 
of four-dimensional simplicial quantum gravity in this article. 

2 Model 

It is still not known how to provide a constructive definition of four- dimensional 
quantum gravity. We have considered a discretized random closed manifold in 
analogy with the two-dimensional case. We numerically evaluated the Euclidean 
path integral with the technique of dynamical triangulation, which gives a discrete 
summing over all possible connections of lattices that may replace the integral 
over diffeomorphism inequivalent metrics. Then, we naturally considered the 
Euclidean Einstein- Hilbert action coupled to Ny copies of U{1) vector fields and 
its discretized model with 4-simplices. The total action is S = Seh + Spi. We 
use the Einstein-Hilbert term for gravity, 

Seh[A, G] = J d'x^ (a - ii?) , (1) 

where A is the cosmological constant and G is Newton's constant. We use the 
discretized action for gravity, 

Seh[k2, Ki] = K4N4 - K2N2 

where K2 ~ l/C; ^4 is related to A' = cA (c is the unit volume) and Ni is the 
number of z-simplices. We use the plaquette action for f/(l) gauge fields |T9| , 

Spi = E [A{h,) + A{l^k) + A{h^)]^ (3) 
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where kj denotes a link between vertices i and j, tijk denotes a triangle with 
vertices i, j and k and o{tijk) denotes the number of 4-simplices sharing triangle 
tijk- We consider the U(l) gauge field A{lij) = —A{lji) on a link lij. 

We consider that a partition function of gravity with A''^ copies of U{1) gauge 
fields is 

Z{k,,^,,Nv)=Y.^"''''' E e'^^^^n/ n dA{l)e-'^'. (4) 

N4 t{2D)£T{4D) Nv 1&{2D) 

We sum over all four-dimensional simplicial triangulation T{4D) in order to carry 
out a path integral over the metric, where we fix the topology with S*^. As is well 
known, we must add a small correction term {6S) |^ to the lattice for fluctuations 
of volume, 

SS = 5k2{N^ - ATf (5) 

where is adjusted with an appropriate choice; we use 5hi2 = 0.00025. 

Near to the critical point, it is expected that the partition function (eq.(H)) 
behaves as 

T(4D) 

~ l/|iV4-A^r''T"^", (6) 

where 'jst is the string susceptibility exponent related to the entropy of the man- 
ifold. In the situation; 7^4 > 0, the manifold grows to the spiky configuration. 
Conversely, in 7^4 < 0, the surfaces grow to smooth structures. 

Using Monte-Carlo simulations, we evaluated the partition function (eq.@). 



We followed the way of updating the configuration after ref. ||T9|. This is that the 
gauge fields renew the sequence of the (p, g)-move according to the weight of the 
Bolzman factor. The action was checked by the Metro police methods. As we held 
the same conditions with Bilke et al., we performed a geometry update, {p,q)- 
move, updating the gauge fields by the heat-bath sweep and the over- relaxation 
sweep. An over-relaxation sweep was introduced for the convergence, 

^ij A.j — '2 Aij. (7) 

The measurement chance (A^4 = iV^"''^'^*) came at intervals of about 100 sweeps, 
where we counted "1 sweep" as the flow of "heat-bath sweep {p, g)-move 
over-relaxation sweep" . The measurement processes and updating processes are 
the same as the pure gravity case. However, since the case of adding vector fields 
costs more CPU time, we arranged for suitable compute time. 



3 Phase diagram with vector fields 

In order to investigate the phase structure of simplicial quantum gravity, first we 
calculated the string susceptibility exponent (7st). This exponent is defined by 
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the asymptotic form of the partition function (eq.(^). It is known that the case of 
7sj corresponds to the dominance of the branched polymer p2| . We measure the 



string susceptibihty exponent (7^4) with the method of the MINBU (Minimum 
Necked Baby Universe) analysis [|I|, 0,|2^. It is a powerful exponent for probing 
the property of quantum geometry, and is easily measured. We can count a baby 
universe which is connected to the mother universe via a minimal neck. The 
distribution function for the MINBU analysis with the size B can written as 

3(A-B + 1)(B + 1)Z[A -B + 1]Z[B + 1] 

= m 

(8) 

where we use the asymptotic form of the partition function, Z[A] ~ A'^^'^'^e'^^, 
where A denotes the volume of the mother universe. 

In the pure gravity measurement of the string susceptibility exponent 




is discussed in ref. From the numerical results, the string susceptibility 

exponent supports the idea that the pure case of simplicial quantum gravity has 
two distinct phases. What is important in the A'^v' = 1 case is that the usual 
phase transition point (/tg) is different from another transition point (^2), which 
separates the •jst < region from the 7^4 > region, and •jst becomes negative 
on the phase transition point (^2). This fact leads to the definition of a new 
smooth phase. This phase is defined by an intermediate region between these two 
transition points and K2). In the pure gravity case, it is clear that 
and thus there is no evidence for the existence of a new smooth phase. On the 
other hand, in the case of Ny = 1 with N4 = 16K, we observe the 7^4 < region 
beyond the usual phase transition point (/tg). We also observe a very obscure 
transition from 7^4 < to 7^4 > at ^2- 

We give the numerical results at Table.|l[ In Fig.||, we plot •jst for various 



Table 1: Measurement results of 7^^ at A^4 = 16K. 
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Nv = 1 


Nv = 2 


1.0 


-0.43(5) 


-2.4(3) 


-3.9(3) 


1.5 


0.46(2) 


-0.87(9) 


-2.4(2) 


2.0 


0.53(1) 


-0.07(3) 


-1.05(9) 


2.5 


0.59(2) 


0.136(4) 


-0.45(5) 


3.0 


0.67(1) 


0.123(3) 


0.01(4) 


3.5 




0.188(7) 


0.15(8) 



numbers of gauge fields versus K2 with volume A^4 = 16i^'. In the case of adding 
vector fields, we can find that the usual phase-transition point (/tg) is different 
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Figure 1: String susceptibility exponents {jst) plotted versus the coupling con- 
stants (^2) for Nv = 0, 1 and 2. 
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Figure 2: Schematic phase diagram of four-dimensional simphcial quantum 
gravity. 



from another transition point (^2) which separates the '^st < region from the 
7sj > region and 7st becomes negative at the phase-transition point (^2)- This 
fact leads to the existence of a new phase. We thus call this intermediate region 
the smooth phase between these two transition points (^2 and n'^. 

We can consider that simplicial quantum gravity coupled to vector fields has 
three phases p0|,|2l|,^. In Fig.H, we show a schematic phase diagram. We have 
three phases in this parameter space: a crumpled phase, a smooth phase (shaded 
portion) and a branched polymer phase. Furthermore, the smooth phase expands 
with adding more vector fields. We show the discontinuous phase transition as 
the thin line and the smooth phase transition as the thick line. In the case of 
adding the vector fields, there are two separate phase transition lines: the usual 
phase transition line {n^ and an obscure phase transition line (^2). The obscure 
transition at ^2 been shown to be third order or cross over which is 



very similar to the c = 1 barrier in two-dimensional quantum gravity [^. In 
two dimensions the c = 1 barrier is well-known as an obscure transition from the 
fractal phase (c < 1 and 7^/^ < 0) to the branched polymer phase (c > 1 and 
ll? > 0)- consider that the obscure phase transition point may be a threshold 
point, like the c = 1 case in two dimensions. From reports of Antoniadis et al. [0, 
the quantum field theory of gravity with conformal invariance has a central charge 
Q^, which has a threshold value, ~ 8. We expect that the obscure phase 
transition point that is defined by the threshold value, 7sf ~ 0, may be the same 
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Figure 3: Node susceptibility (x) plotted versus the coupling constant {1^2) for 
Nv = 1. 



situation as the ~ 8 case. Another phase transition point is the usual phase 
transition point between the strong coupling region (the crumpled phase) and the 
intermediate region (the smooth phase). We expect that the phase transition at 
K2 is continuous. This leads to the continuum limit of four- dimensional quantum 
gravity. 

Now, let us take a look at the transition at K2 in the case of Ny = 1. In 
order to investigate the phase transition, we observe at the exponents of the 



node susceptibility. The node susceptibility is defined in ref. |]T2[ as follows: 



X = ^ (< iVo^ > - < iVo >') . (9) 

In Fig.|^, we plot the node susceptibihty(x) as a function of K2 with the volume 
A^4 = 16K, 2AK and 2>2K, respectively. The node susceptibility (x) has a peak 
value at the critical point (^2). We find the peak value in each size. The height 
and the width of the susceptibility peaks give the finite size scaling exponents of 
the phase transition. The peak value {xmax) and the width of peak (5^2) grow 
as Ni in power. The susceptibility exponents (A and V) are defined by fTT 



XmaxOCiVf, {6k2(xN^'). (10) 
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Figure 4: Histogram of node (A^o) near to the critical point {k,2 = 1.37147(1)) 
for Ny = 1 with A^4 = 32K. The number of configurations is 10,000. 



From the numerical results (Fig.^, we obtain the susceptibility exponents: A = 
0.4(1), r ~ 0.5(3). These values are apparently smaller than 1, though they 
are 1 |12[ in the pure gravity case. This numerical results show that simplicial 



quantum gravity coupled to a vector field has the different type of phase transition 
from simplicial quantum gravity in the pure gravity case. 

In Fig.§, we show a histogram of Nq for a size of N4 = 32K near to the 
critical point (kj = 1.37147(1)). In the pure gravity case, previously, a double 
peak structure has been found [jlT],|l2[. This structure relates to the latent energy. 
The fact that the phase transition is first order is shown. However, by adding 
the vector field, the double peak structure disappears. We thus conclude that 
the phase transition between the crumpled phase and the smooth phase may be 
continuous, not first order. 

4 Scaling property of four-dimensional 
simplicial quantum gravity 

Recent numerical results obtained by dynamical triangulation in four space-time 
dimensions suggest the existence of the scaling behavior, for example, the MINBU 



distribution ||2^ , correlation functions ||2^ and the boundary volume distributions 



rO[] . If we assume the existence of a correct continuum limit, the scaling property 
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of quantum geometry gives important information about the continuum theory, 
because some scahng properties are related to the universahty of the theory. One 
of the interesting observes is the fractal dimension (Hausdorff dimension) (dn), 

This is based on studying the behavior of the volume V^^\D), the number of 
simplices, within a geodesic distance D. The geodesic distance D is defined as 
the shortest distance between two simplices through the center of the simplices. 
In the pure gravity case, the Hausdorff dimension shows a different behavior in 
each phase. In the crumpled phase, the Hausdorff dimension diverges {dn — >■ oo). 
The Hausdorff dimension rises very steeply to a large value below the transition, 
while it falls very rapidly to a branched polymer, dn = 2, above the transition. 
However, in the case of gravity coupled to vector fields, each of three phases 
shows a different behavior for the Hausdorff dimension. In the crumpled phase 
and the branched polymer phase, the behavior for the Hausdorff dimension is 
similar to that of pure gravity. Furthermore, we find that the smooth phase is 
the intermediate region with 2 < rfji/ ~ 4. For the smooth phase, the value of 
the Hausdorff dimension (dn) is changed to smooth, as compared with the pure 
gravity case. This fact supports the results of a finite size analysis. Especially, 
at the critical point, we observe the Hausdorff dimension, dn = 4.6(2), with 
N4 = 32K. 

Next, let us discuss the scaling structure of four- dimensional DT mfd, at the 
focusing on the scaling structure of the boundaries in four-dimensional Euclidean 
space-time using the concept of geodesic distances. We consider that the scahng 
structure of the boundaries has more informations than the Hausdorff dimension, 
and that it is the way of directly searching for the structure of quantum geometry. 
In the pure gravity case, the scaling property of the boundary volume is discussed 
in ref . fl^ . As in the previous analysis 0] , we assume that the boundary volume 



distribution {p{x, D)) is a function of a scaling variable in analogy with the loop 
length distribution function in the two-dimensional model. Fortunately, in the 
two-dimensional model, the loop length distribution function has been calculated 
analytically The loop length distribution function, p{x = L/D^,D), which 
gives the probability of the boundary loop with the loop length (L) within a 
geodesic distance (D), is given as a function of the scaling variable x = L/D"^: 

p{x = LID\D) = + + Y^'^'Y'^- (12) 

This distribution function consists of two different types of distributions. The 
first two term of eq.(0) represent the baby loops that the universe has a small 
boundary volume; the last term represents the mother loop that the universe has 



a large boundary volume. The distribution function, p(x = L/D ,D) of eq.(|T^ 
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Figure 5: Schematic picture of the boundaries (Si,S2 and S3) at distance D 
and a boundary (Eg) at distance D' in four- dimensional Euchdean space with 
topology. The mother boundary is defined as a boundary (E2) with the largest 
tip volume (V2); the other boundaries are defined as baby boundaries. 



satisfies the scaling relation under rescaling (D D' = ^/XD, L ^ L' = XD): 

p{L,D) = X-'p{L',D'). (13) 

In the four-dimensional model, unfortunately, a similar scaling relation about 
the boundary volume is not yet known. We thus assume the distribution function, 
p{x = V/D"'), of the volume of the boundary {V) within the geodesic distance 
[D) in four- dimensional dynamical triangulated manifold. This is a function of 
the scaling variable, x = V/D"", with scaling parameter [a) in analogy with the 
two-dimensional model. In Fig.^ we show a schematic picture of our boundary 
analysis. 

First, we consider the scaling structures of these three phases: a crumpled 
phase, a smooth phase and a branched polymer phase ||20|. Actually, in the 



smooth phase we observe that the distribution (p) becomes fractal in the sense 
that the sections of the manifold at different distances from a given 4-simplex look 
exactly the same after a proper rescaling of the boundary volume. Furthermore, 
the shape of this scaling function is very similar to that of the two-dimensional 
case 0,0. The best account for this excellent agreement in the four-dimensional 
case can be found in the dominance of a conformal mode and a fractal property. 
We have also investigated the boundary volume distribution in the crumpled 
phase and the branched polymer phase. It seems reasonable to suppose that this 
new smooth phase has a similar fractal structure to that of the two-dimensional 



10 



fractal surface, and that there is a possibihty of taking a continuum hmit in the 
phase. We have also investigated the boundary volume distribution in both the 
crumpled phase and the branched polymer phase. In the crumpled phase we find 
that one mother universe is dominant, while in the branched polymer phase there 
is no evidence for the existence of a mother universe. 

Next, let us discuss the relation between the scaling parameter (a) and 
the Hausdorff dimension {(!}{)■ The expectation value of the boundary three- 



dimensional volume appearing at distance D has been introduced in ref. [10 



1 r'^ 

< -1/(3) >= — / dV Vp{x = V/D", D), (14) 

where Vq denotes the UV cut-off of the boundary volume and is the normaliza- 
tion factor. If the boundary volume has the scaling property with the universal 
distribution {p{x,D)) and Vq — > 0, 

< >^ (15) 

Then, we should obtain a finite fractal dimension, 

df = a + l, (16) 

with the fractal dimension df. We measure the volume of the mother boundary 
as a function of D. The mother boundary is defined by the boundary having the 
largest tip volume. In Fig.^, we plot the mother boundary volume < V^^^ > with 
the size of A^4 = 32K at the critical point. As a result, the mother boundary 
volume shows a scaling, and we obtain the scaling parameter (a = 3.7(5)). Then, 
we can estimate the fractal dimension {df = 4.7(5)). On the other hand, we 
measured the Hausdorff dimension, which results in dn = 4.6(2). Both results 
are consistent {df ^ d^)- We thus expect that the boundary volume has a 
scaling property in the sense that the manifold at different distances from a given 
4-simplex looks exactly the same after a proper rescaling of the boundary volume. 



5 Modified balls-in-boxes model 

In this section, we discuss a scenario for the phase structure and the phase tran- 
sition in four-dimensional simplicial quantum gravity coupled to matter fields. In 
the pure gravity case, the first order phase transition is described by the balls-in- 
boxes model p6|-p!8|. This model is considered to be a simple mean field model 
about simplicial quantum gravity. It describes a fixed number balls distributed 
into a variable M of boxes. The partition function is given by 

Zm,N = PiQl) ■ ■ ■Pi9M)^qi+-+<lM,N, (17) 
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Figure 6: Power fit to the mother boundary three-dimensional volume < V^^"^ > 
at the critical point (k^ = 1.37147(1)) for Ny = I with = 32K. 

where is the number of balls in box i, and p{q) denotes a weight which depends 
on only the number of balls (q) in the box. The original balls-in-boxes model has 
the following correspondence: 

Vertex ^ Box and Simplex ^ Ball. 

Therefore, qi represents the vertex order. The original model has been analyzed 
in the case of power like weights: p{q) = q~^ (when /5 is a parameter for a weight). 
The system has two distinct phases: a crumpled phase and an elongated phase, 
and it has a discontinuous phase transition for (3 > 2. It describes the phase 
structure about four- dimensional pure simplicial quantum gravity. However, this 



model describes the continuous phase transition at the (5 < 2 ||18[. We have 
noticed this fact, which gives one of the motivations for the modified model of 
simplicial quantum gravity. In the model for simplicial quantum gravity, we can 
consider that parameter j3 is related to the measure term. If we consider the 
conformal matter fields (as the conformal charge), parameter j3 has the same 
origin as the conformal anomaly, because they are related to the measure term. 



The gauge field gives rise to a modified measure factor |2T 

Mt = \{o{tf. (18) 

t 

This model shows three distinct phases: a crumpled phase, a crinkled phase 
(correspond to a smooth phase) and a branched polymer phase. It is consistent 
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with our results. The natural correspondence between our phase structure and 
that in Ref. is as follows: the transition between the crumpled phase and the 



branched polymer phase is discontinuous (probably, first-order phase transition), 
the transition between the crumpled phase and the crinkled phase is continuous, 
and the transition between the crinkled phase and the branched polymer phase 
is obscure. 

Thus, in order to investigate the phase structure of adding matter fields, we 
exchange the relation of "vertices-simplices" into that of "triangle-simplices" : 

Triangle Box and Simplex <->■ Ball. 

We thus introduce the triangle order (o(t)) instead of the vertex order (g). This 
modification causes the partition function to give rise a constant shift for the 
constraint: 

E = 5iV4, (19) 

i:vertex 

E o{t) = lONi. (20) 

t-.triangle 

We analyzed the modified balls-in-boxes model with using the same method of 
the both original models. 

We thus considered the case of simplicial quantum gravity coupled to matter 
fields. The partition function is 

P ; -^matter ) ^ ^ ^ ' 

m 

E P(l^) ■ ■ ■p(gm)5gi + ...+g„,iv(^ma«ar(m))^'— ^(21) 

Now we discuss the following two cases: (a) the case of coupling to Nb copies of 
the boson fields and (b) the case of coupling to Ny copies of the vector fields, (a) 
In the first case, the matter action is 

Sboson = ^{(j>i - <f)jY. (22) 
(ij) 

We discuss the effect from the matter. From the perturbation, the Gaussian 
boson matter gives 

Zboson = E Constant. (23) 

triangles 

We can estimate the effect to the weight of one-boxes: 

p{q)^C-p{q). (24) 

The Gaussian bosons give the effect of multiplication by a constant. From a mean 
fields analysis, the phase diagram will be changed. However, it is known that the 
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phase diagram is not changed by only a few Gaussian boson fields, according 
to numerical analysis [^. We consider that the multiplication by a constant is 



not enough to change the phase diagram for "real" simplicial quantum gravity, 
because of the too small effect of multiplication by a constant, (b) In the second 
case, we consider adding the Gaussian U{1) gauge matter. If we use the plaquette 
action for the gauge fields, the partition function is written as 

Zgauae= E /(oW), (25) 

triangle 

where / is a function of the triangle order (o(t)). From the perturbation about 
the gauge interaction, f{o{t)) = Co + Cio{t). Then, the partition function is 
replaced, as follows: 

Zgauge ^^0{t). (26) 

t 

This is just the modified measure. We estimate the effect from the gauge matter 
fields as 

(27) 

That is to say, the Gaussian gauge matter fields lower f3. 



6 Summary and Discussion 

We have investigated the phase structure and the phase transition with a model 
of four-dimensional simplicial quantum gravity coupled to U{1) gauge fields. 

The results of our study are summarized by the schematic phase diagram in 
Fig.|^. We checked this phase diagram in the case of Ny = 1, 2, 3 at a volume of 
A^4 = 16K. We found three phases in this parameter space: a crumpled phase 
(7sj —^oo,df^ oo), a smooth phase (7^* < 0, 2 < ~ 4, N0/N4 < 0.25) and a 
branched polymer phase (7^4 > 0, df = 2, Nq/N^ = 0.25). The thin line denotes 
a discontinuous phase-transition line which is known in pure gravity; moreover, 
the thick line denotes a smooth phase-transition line. As contrasted with the 
phase diagram of pure gravity, the phase diagram means richer structures. In the 
crumpled phase one can find singular sub-simplices, for example, vertex order 
concentration and link order concentration. The smooth phase is defined as a 
region between the critical point (^2) and the obscure phase transition point 
(K2). We observe the negative string susceptibility in this region. Then, with 

= IQK {k2 = 1.7) and Ny = 1, we obtained 7,* = -0.38(5), dj = 2.8(5) and 
a good scaling relation of the boundary volume distributions. We consider that 
the scaling structure of this smooth phase is similar to that of a two-dimensional 
random (fractal) surface. This smooth phase is slowly broken to a branched 
polymer phase that has no mother structure at the obscure phase transition 
(^2). We obtained an obscure-transition line (a broken line in Fig.|^); therefore. 
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we suggest that the obscure-transition corresponds to the c = 1 barrier in two- 
dimensional quantum gravity. 

As for the phase transition at the critical point (k2 = K2), we showed the finite 
size scaling at the critical point and the histogram of the node. We calculated 
some finite size scaling exponents, and showed that the value of these exponents 
is smaller than 1. It has been discussed [0 that the value of 1 is expected to 
be a first order phase transition for the results of pure gravity. However, in the 
case of adding one gauge matter field {Ny = 1), the numerical results show that 
the phase transition is smooth, in the contrast to pure gravity, and then that a 
double peak structure of the node histogram disappears. 

Furthermore, in order to investigate the property of quantum geometry and to 
discuss the universality of the manifold, we considered the scaling property of the 
boundary volume. In the smooth phase with A^4 = 16K {k2 = 1.7) and Ny = 1, 
we obtained df = 2.8(5) and a good scaling relation of the boundary volume 
distributions. The scaling structure of this smooth phase is similar to that of a 
two-dimensional random (fractal) surface. This suggests the existence of a new 
smooth phase in four-dimensional simplicial gravity. The other two phases have 
a similar scaling property to that of pure gravity. We have shown the scaling 
property of the mother boundary, where the scaling parameter is consistent with 
the Hausdorff dimension. Then, the boundary volume has a scaling property 
with the scaling variable (x = V/D'^f^'^). We expect that the boundaries have a 
fractal structure and universality of the scaling relations. We also discussed the 
modification of the balls-in-boxes model. The role of a vertex is exchanged into a 
triangle. This clarifies the relation between the measure factor of the numerical 
model and that of the analytical mean field model. 

We expect the existence of genuine four-dimensional quantum gravity on the 
critical point (/tj) with the vector fields. Our recent investigations will give further 
evidence for the existence of an ultraviolet fixed point of the quantum field theory 
of gravity. 
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